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ABSTRACT 

We report the discovery of two super-Earth mass planets orbiting the nearby K0.5 dwarf HD 7924 
which was previously known to host one small planet. The new companions have masses of 7.9 and 
6.4 M 0 , and orbital periods of 15.3 and 24.5 days. We perform a joint analysis of high-precision radial 
velocity data from Keck/HIRES and the new Automated Planet Einder Telescope (APE) to robustly 
detect three total planets in the system. We refine the ephemeris of the previously known planet 
using five years of new Keck data and high-cadence observations over the last 1.3 years with the APE. 
With this new ephemeris, we show that a previous transit search for the inner-most planet would have 
covered 70% of the predicted ingress or egress times. Photometric data collected over the last eight 
years using the Automated Photometric Telescope shows no evidence for transits of any of the planets, 
which would be detectable if the planets transit and their compositions are hydrogen-dominated. We 
detect a long-period signal that we interpret as the stellar magnetic activity cycle since it is strongly 
correlated with the Ga H H & K activity index. We also detect two additional short-period signals 
that we attribute to rotationally-modulated starspots and a one month alias. The high-cadence APE 
data help to distinguish between the true orbital periods and aliases caused by the window function 
of the Keck data. The planets orbiting HD 7924 are a local example of the compact, multi-planet 
systems that the Kepler Mission found in great abundance. 

Subject headings: planetary systems — stars: individual (HD 7924) 


1. INTRODUCTION 

The archetypal planets of our Solar System—Jupiter 
the gas giant, Neptune the “ice” giant, and Earth the ter¬ 
restrial planet—represent an incomplete inventory of the 
planet types in our galaxy. We are locally impoverished 
in “super-Earths,” the broad category of planets inter¬ 
mediate in size and mass between Earth and Neptune. 
Doppler searches of nearby stars showed tha t super- 
Earth planets in close-in orb its are plentiful (Howard 
etll][m 0 |Mayor et al.||2011|). Results from the Kepler 
mission confirmed and rehned our knowledge of the size 
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(Howard et al. 2013 Pep e et al. 2Q13|) to pu ffy planets 
with thick gas envelopes ( [Kipping et al.||2014 ). Interme¬ 
diate planets with densities of ^ 3 g cm“^ are consistent 
with a broad range of interior structures and atmosphere 
sizes. Planets smaller than ^1.6 Earth radii (Re) are 
more likely t o have a high density and presume d rocky 
composition ( Weiss fc Marcy|[2Q14 Rogers||2Q15 ). 

The large population of super-Earths orbiting close 
to their host stars was a surprise. Population synthe¬ 
sis models of planet fo rmation had predicted that such 
syste ms would be rare ( Ida fc Lhi]|2QQ4 Mordasini et al. 


2009). Planet cores were expected to mostly form be¬ 
yond the ice line and rarely migrate to close orbits un¬ 
less they first grew to become gas giants. Nevertheless, 
close-in, low-mass planets are com mon and often appear 


in compact multi-pla net systems (Lissauer et al. 2011 


Pang & Margot 2012). Theoretical models are catching 


up, with refinements to the disk migration and multi- 


models (|Ida & Lin||2010: 

Alibert et al.l 

2013 

Schlichting 

2014 

Lee et al. 2014|). i 

\ new class ot “in situ” torma- 


tion niodels have also been proposed in which systems of 
super-Eart hs and Neptnnes emerge naturally from mas¬ 


sive d isks (Hansen & Murray 2012 Ghiang & Laughlin 


2013). 


The Eta-Earth Survey (Howard et al. 2010) at Keck 
Observatory played a important role in the discovery that 
super-Earths are abundant. Using the HIRES spectrom¬ 
eter, our team searched for planets in a volume-limited 
sample of 166 nearb y G and K dwarf stars. Our search 
yielded new planets ( Howard et al.||200'9 2011a|b 2014) 
and detection limits for each star. Putting these to¬ 
gether, we measured the prevalence of planets in close-in 
orbits as a function of planet mass (Msin/). This mass 
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function rises steeply with decreasing mass: planets in 
the mass range 3-10 are about twice as common as 
10-30 M 0 planets. 

The first low-mass pla net discovered in the Eta-Earth 
Survey was HD 7924b (Howard et al. 2009, H09 here¬ 
after), a super-Earth with a mass of 8.7 M 0 and an 
orbital period of 5.4 days. We have continued to ob¬ 
serve HD 7924 and other Eta-Earth Survey stars using 
HIRES. These additional measurements probe smaller 
masses and larger star-planet separations. We have also 
started observing a subset of the Eta-Earth Survey stars 
with t he Automated Planet Einder (APE; |Vogt et al. 
2014a), a new telescope at Lick Observatory. APh' is a 
robotic 2.4-m telescope designed exclusively for Doppler 
discovery of ex oplanets. It feeds the hi gh-resolution Levy 
Spectrometer (Radovan et al. 2010) that uses an io¬ 
dine reference spectrum to c alibrate tne wavelen gth scale 
and point spread function (Butler et al. 1996), achiev¬ 
ing a Doppler precision similar to HIRES while running 
without human intervention during an observing night. 
APE exploits high measurement cadence (nearly nightly) 
to disentangle the complicated low-amplitude signals of 
multi-planet systems in the face of stellar activity. 

In this paper we announce two additional super-Earths 
orbiting HD 7924 based on RVs from the Keck-HIRES 
and APE-Levy spectrometers. We describe the proper¬ 
ties of the star HD 7924 in Sec. and our Doppler mea¬ 
surements from APE/Levy and Keck/HIRES in Sec. 1^ 
Our analysis of the RV data, including discovery of the 
three Keplerian signals and consideration of false alarm 
probabilities, alias periods, and chromospheric activity, 
is described in Sec.lU We conclude with a discussion and 
summary in Sec. 

2. STELLAR PROPERTIES 

HD 7924, also known as HIP 637 9 or GJ 56.5, is a 
nearby (16. 82 pc; van LeeuwenI 200 7) and bright K0.5V 
dwarf star (von Braun et aL||2(J14p. It is slightly metal 
poor relative to the Sun and nosts one previously known 
planet. It is relatively in active with logR^K = —4.89 
(Isaacson & Eischer 2010), but we do detect some low- 
leve l chromospheric activity nonetheless (see Section 
T.5|). We list our adopted stellar parameters in Table 


Most of our spectroscopically-derived stellar parame¬ 
ters are adopted from HQ9, wh ich were originally de¬ 
rived by Valenti & Eischer (2005) using the Spectroscopy 
Made Easy (SME) LT'E spectral synthesis code. How¬ 
ever, HD 7924 has been the fo cus of several studies since 
the discovery of HD 7924b. Santos et al. (2013) per¬ 
formed a uniform analysis or 48 planet-hosting stars. 
They find Teff = 5133 ± 68 K, \ogg = 4.46 T 0.12. 


and [Fe/H] = —0.22 ± 0.04 for this st ar, von Braun 


et al. (2014) use the empirical relations of Boyajian et al 


( |2012 ) to calculate a mass of 0.81 M 0 for HD 7924 with 


an error estimate of 30%. We computed the stellar mass 


and radius from Teff, log^f, and [Fe/H] using the Tor 
res et al.|(| 2010 |) relations and found 0.81 ± 0.02 M 0 and 
0.75 ± 0.03 R© . All of these values are within 1-a of our 
adopted values. 


HD 7924 was observed by von Braun et al. (2014) 
using long-baseline interferometry on the Georgia State 
University Genter for High Angular Resolution As tron- 
omy (CHARA) Array ( ten Brummelaar et al.|20Q5 ). We 


TABLE 1 

Adopted Stellar Properties of HD 7924 


Parameter HD 7924 


Source 


Spectral type 
B — V (mag) 

V (mag) 

J (mag) 

H (mag) 

K (mag) 
Distance (pc) 
Teff (K) 
log^f (cgs) 
[Fe/H] (dex) 
vsini (kms“^) 

L. (To) 

M, {Mq) 

R, (To) 

log 

8'hk 


K0.5V 

0.826 

7.185 

5.618 ± 0.026 
5.231 ± 0.033 
5.159 ± 0.020 
16.82 ± 0.13 
5075 ± 83 


4.59 


+0.02 

- 0.03 


-0.15 ± 0.03 
1.35 ±0.5 
0.3648 ± 0.0077 


0.832 


+0.022 

- 0.036 


0.7821 ± 0.0258 


-4.89 


0.20 


Ivon Braun et al.l (|2014| 


TTug- 

H09 


Cutri et al. 

(2003 


Outri et al. 

(2003 


Cutri et al. 
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Takeda et al. ( 20071 
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2014) 

Isaacson & L'lsclier 

(2016 

Isaacson & L ischer 
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adopt their value for Teff of 5075 ± 83 K that they ob¬ 
tained by fitting the spectral energy distribution from op¬ 
tical through infrared wavelengths. They also obtained 
precise values for the luminosity and radius that we list in 
Tablelil These values are consistent within 1 — a with the 
values adopted in H09 that were based on fitting stellar 
evolut i on mo dels with spectroscopic parameters. |M ason| 


et al. ( 2011 ) conducted high-resolution speckle imaging 
of HD 7924 and detected no companions within three V 
magnitudes of HD 7924 with separations between 0703 
and 175. 


3. MEASUREMENTS 
3.1. Keck/HIRES Spectroscopy 

We collected 599 new RV measurements for HD 7924 
over the last 5 years since the discovery of HD 7924b 
(H09) using the HIRES spectrograph on the Keck I tele¬ 
scope (|Vogt et al. 1994). When this new data is combined 
with the data from H09 and APE we have over 10 years 
of observational baseline (see Figure!^). Our data col¬ 
lection and reduction techniques are aescribed in detail 
in H09. On Keck/HIRES we observe the star through a 
cell of gaseous iodine in order to simultaneously forward 
model the instrumental line broadening function (PSF) 
and the subtle shifts of the stellar lines relative to the 
forest of iodine lines. The HIRES detector was upgraded 
in August of 2004 and the RV zero-point between the 
pre-and post-upgrade data may not necessarily be the 
same. For this reason, we allow for separate RV zero- 
points for the pre-upgrade, post-upgrade, and APE data 
sets. All RV measurements and associated Ga II H & K 
5'hk activity indices (the ratio of the flux in the cores 
of the H & K lines to neighboring continuum levels) are 
listed in Tabled 


3.2. APF/Lcvy Spectroscopy 

The Automated Planet Finder is a 2.4 m f/15 
Cassegrain telescope built by Electro-Optical Systems 
Technologies housed in an IceStorm -2 dome located at 
Lick Observatory atop Mount Hamilton, 20 miles east of 
San Jose, California. The telescope operates completely 
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unattended using a collection of Python, Tel, bash, and 
csh scripts that interact with the lower level softwa re op¬ 


erating on the Keck Task Library keyword system (Lup- 


ton & Conrad 1993). Every scheduled observing night. 


our automation software queries an online Google spread¬ 
sheet (that serves as our target database) and creates 
two observing plans: one for good conditions and the 
other for poor. Weather permitting, shortly after sun¬ 
set the observatory opens automatically and determines 
whether to use the observing plan designed for good or 
poor conditions by monitoring a bright star. Once the 
observing plan is started, the high-level software contin¬ 
ues to monitor the conditions and adjusts the schedule 
or switches observing plans if necessary. Observations of 
our targets (primarily for radial velocity measurements) 
continue throughout the night until conditions deterio¬ 
rate too much to remain open or the until morning 9- 
degree twilight. 

The Levy Spectrograph is a high-resolution slit-fed op¬ 
tical echelle spectrograph mounted at one of the two 
Nasmyth foci of the APE designed specific ally for the 


detect i on and characterization of exoplanets (Vogt et al. 
2Q14b; Radovan et al.|2Q14 Burt et al.]2Q14 ). Eadi spec- 


trum covers a continuous wavelength range from 3740 A 
to 9700 A. We observed HD 7924 using a I'.'O wide decker 
for an approximate spectral resolution of R = 100, 000. 
Starlight passes through a cell of gaseous iodine that 
serves as a simultaneous calibration source for the in¬ 
strumental point spread function (PSE) and wavelength 
reference. In addition, we collected a high signal-to-noise 
spectrum through the 075 wide decker {R = 150,000) 
with the iodine cell out of the light path. This spectrum 
serves as a template from which we measure the rela¬ 
tive doppler shifts of the stellar absorption lines with 
respect to the iodine lines while simultaneously mod¬ 
eling the PSE and wavelength scale of each spectrum. 
APE guides on an image of the star before the slit using 
a slanted, uncoated glass plate that deflects 4% of the 
light to a guide camera. Phot on-weighted times of mid¬ 
exposure are recorded using a software-based exposure 
meter based on guide camera images that monitors the 
sky-subtracted light entering the slit during an exposure 
( Kibrick et al.||2QQ6 ). 

We measure relative radial velocities (RVs) using a 
Doppler pipeline descended from the iodine technique in 
Butler et al. (1996). Eor the APE, we forward-model 848 


segments of each spectrum between 5000-6200 A. The 
model consists of a stellar template spectrum, an ultra 
high-resolution Eourier transform spectrum of the iodine 
absorption of the Levy cell, a spatially-variable PSE, a 
wavelength solution, and RV. We estimate RV uncertain¬ 
ties (Table[^ as the uncertainty on the mean RV from the 
large number of spectral segments. Well-established RV 
standard stars have a RV scatter of ^2-3 ms“^ based on 
APE measurements. This long-term scatter represents a 
combination of photo-limited uncertainties, stellar jitter, 
and instrumental systematics. We collected a total of 
109 RV measurements of HD 7924 on 80 separate nights 
(post outlier rejection) over a baseline of 1.3 years with 
a typical per-measurement uncertainty of 2.0 m s“^. 

Eor both the APE and Keck data, Julian dates of the 
photon-weighted exposure mid-times were recored dur¬ 
ing the observations, then later converted to Barycentric 


TABLE 2 

Radial Velocities of HD 7924 


BJDtdb 
(- 2440000) 

RV 

(ms-i) 

Uncertainty 

(ms-l) 

Instrument^ 

8'hk 

I2307.77I62 

5.23 

I.3I 

k 


12535.95639 

3.24 

1.16 

k 


13239.08220 

-3.00 

0.91 

j 

0.218 

13338.79766 

3.83 

1.08 

j 

0.227 

I6505.9973I 

-7.97 

1.46 

a 


I65I5.90636 

3.88 

1.63 

a 



(This table is available in its entirety in a machine-readable form in 
the online journal. A portion is shown here for guidance regarding 
its form and content.) 

^ k = pre-upgrade Keck/HIRES, j = post-upgrade Keck/HIRES, 
a = APE 


Julian date in t he dynamical time sys tem (BJDtdb) us¬ 
ing the tools of Eastman et al. I ( |2010[ P 


We rejected a handful of low signalTo-noise ratio spec¬ 
tra (S/N < 70 per pixel) and measurements with un¬ 
certainties greater than nine times the median absolute 
deviation of all measurement uncertainties relative to the 
median uncertainty for each instrument. This removed 
a total of 11 RVs out of the 906 total measurements in 
the combined data set (Keck plus APE). 

After outlier rejection we bin together any velocities 
taken less than 0.5 days apart on a single telescope. 
Since data taken in short succession are likely affected by 
the same systematic errors (e.g. spectrograph defocus), 
these measurements are not truly independent. Binning 
helps to reduce the effects of time-correlated noise by 
preventing multiple measurements plagued by the same 
systematic errors from being given too much weight in 
the Keplerian analysis. When the data are binned to¬ 
gether and the uncertainty for the resulting data point is 
added in quadrature with the stellar jitter (an additional 
error term that accounts for both stellar and instrumen¬ 
tal systematic noise), the binned data point receives only 
as much weight as a single measurement. While this 
likely reduces our sensitivity slightly we accept this as a 
tradeoff for more well-behaved errors and smoother 
surfaces. An i ndep endent analysis using the unbinned 
data in Section |4.4| finds three planets having the same 
orbital periods, eccentricities, and masses within 1-cr as 
those discovered by analyzing the binned data. 


4. KEPLERIAN ANALYSIS 
4.1. Discovery 

We identify significant periodic signals in the RVs using 
an iterative multi-planet detection algorithm based on 
the two-dime nsional Keplerian Lo mb-Scargle (2DKLS) 
periodogram (O’Toole et al. 2009). Instead of fitting si¬ 
nusoidal func tions to the Rv tirne series (|Lomb||1976t 
Scargle 1982), we create a periodogram by fitting the 
RV data with Keplerian orbits at many different starting 
points on a 2D grid over orbital period and eccentricity. 
This technique allows for relative offsets and uncertain¬ 
ties between different data sets to be incorporated di¬ 
rectly into the periodogram and enhances the sensitivity 


IDL tools for time systems conversion; 
http://astroutils.astronomy.ohio-state.edu/time/ 
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to moderate and high eccentricity planets. 

We fit the Keplerian models using the Levenberg- 
Marquardt (L- M) y^-minimization routi ne in the RVLIN 
IDL package (jWright & Howard] |2009|). Multi-planet 
models are sums ot single-planet models, with planet- 
planet gravitational interactions neglected. Such an ap¬ 
proximation is valid since the interaction terms are ex¬ 
pected to be <Cl ms“^for non-resonant, small planets. 
We def ine a grid of search period s following the prescrip¬ 
tion of Horne & Baliunas (1986) and at each period we 
seed an L-M ht at hve evenly-spaced eccentricity values 
between 0.05 and 0.7. Period and eccentricity are con¬ 
strained to intervals that allow them to vary only half 
the distance to adjacent search periods and eccentricities. 
All other model parameters are free to vary, including the 
parameters of any previously identified planets. The pe¬ 
riod and eccentricity for previously identified planets are 
constrained to be within ±5% and lio%: respectively, 
of their initial values but all other parameters are un¬ 
constrained. This prevents slightly incorrect fits of the 
first detected planets from injecting periodic residuals 
that could mimic further planetary signals. The 2DKLS 
periodogram power at each point in the grid is 


Z(P,e) = ^" 

Xb 


( 1 ) 


where is the sum of the squared residuals to the cur¬ 
rent + 1 planet fit, and x% is the sum of the squared 
residuals to the best Wplanet fit. In the first iteration 
of the planet search (comparing a 1 -planet model to a 
0 -planet model), x% is the sum of the squared error- 
normalized residuals to the mean (the B subscript stands 
for baseline). The 2D periodogram is collapsed into a 
ID periodogram as a function of period by taking the 
maximal Z for each period searched (i.e. the best fit 
eccentricity for every period). 

We start the iterative planet search by comparing a 
0 -planet model (flat line) to a grid of 1 -planet models. 
A strong signal at a period of 5.4 days is detected at 
very high significance in the first iteration. This planet 
was initially discovered by H09 using only 22 % of the 
data presented in this work. When we calculate two vs. 
one planet and three vs. two planet periodograms, we 
find two additional highly significant signals with peri¬ 
ods of 15.3 days and 24.5 days respectively (see Figure 
[^. These two periodic signals are best fit by Keplerian 
orbital models with semi-amplitudes of 2.3 ms“^ and 1.7 
ms“^, respectively, with no significant eccentricity. See 
Tables m and [U for the full orbital solution. 

We find a fourth, significant signal at ^2400 days, but 
we also find a very similar signal (approximately sinu¬ 
soidal with the same period and phase) upon inspection 
of the time series of 5 'hk stellar activity measurements 
(see Figure Although we fit for this periodic sig¬ 
nal as an adaitional Keplerian we do not interpret this 
as an additional planet. Instead we interpret this as 
the signature of the stellar magnetic activity cycle. We 
searched for a fifth periodic signal and found two addi¬ 
tional marginally significant peaks at 17.1 days and 40.8 
days. Due to their marginal strengths and the fact that 
the two periods are related by the synodic month alias 
(1/17.1 ^ 1/40.8 + 1/29.5) we are especially cautious in 
the interpretation of these signals. We scrutinize these 


TABLE 3 

Orbital Parameters 


Parameter 

Value 

Units 

Modified DE-MCMC Step Parameters^ 

log(Pt,) 

0.73223 ±2e - 05 

log(days) 

Vet cos 0)6 

0 15 +0-13 
-0.17 


Ve6sin 0)6 

-0 09 

u.un _o.i5 


log(A'6) 

0.555 

ms“^ 

log(Fc) 

1.184663 

log(days) 

•y/ej cos O)o 

0 20 

U.ZU _Q 24 


y/ejsino)o 

-0.20 


\og(Kc) 

0 364 

U.dOT _0.040 

m 

log(Pd) 

1 QSSQ +0.00027 
i.OOOO _o. 00031 

log(days) 

V^COSO)d 

f) Q1 +0.16 
-0.24 


V^smo)d 

0 02 

U.UZ _Q 37 


\og{Ka) 

0 219 

-0.057 

m 

Model Parameters 


Pb 

5.39792 ±0.00025 

days 

T:onj ,£» 

2455586.38 lo.no 

BJDtdb 

eb 

n +0.056 

U.055 _0.040 


UJb 

339 

-50 

degrees 

Kb 

3 59 

o.on _o.i9 

m s“^ 

Pc 

15 299 +0 0032 
-0.0033 

days 

T 

conj ,c 

2455586.29 

BJDtdb 

ec 

0.098 


LOc 

97 +52 
-60 

degrees 

Kc 

9 +0-21 
-0.20 

m s“^ 

Pd 

24 451 +0-015 
zn.noi _o.oi 7 

days 

Pconijd 

2455579.1 11° 

BJDtdb 

e-d 

0 21 +0-13 
-0.12 


^d 

119 lg7° 

degrees 

Kd 

1.65 ±0.21 

m s“^ 

Tpost-upgrade Keck 

-0.19 ±0.16 

m s“^ 

Tpre-upgrade Keck 

2.0 Ri 

m s“^ 

7APF 

0.28 1°;R 

m s“^ 

c^jitt 

2-41 IRJ 

m s“^ 


^ MCMC jump parameters that were modified from 
the physical parameters in order to speed conver¬ 
gence and avoid biasing parameters that must phys¬ 
ically be finite and positive. 

candidate signals in depth in Sections |4.4| and |4.5[ We 
also see a peak at around 40 days in a periodogram of 
the 5'hk time series after the long-period signal from the 
stellar magnetic activity cycle is removed that is presum¬ 
ably the signature of rotationally modulated star spots. 
We conclude that the 17.1 and 40.8 day signals are most 
likely caused by rotational modulation of starspots, but 
the 5.4, 15.3, and 24.5 day signals are caused by three 
planetary companions with minimum (M^sinip) masses 
of 8.7 M 0 , 7.9 M 0 , and 6.4 M^. 

4.2. Characterization 

We determine the orbital parameters and associated 
uncertainties for the three planet system using the ExoPy 
Differential-Evolution Markov Chain Monte Carlo (DE- 
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Fig. 1.— Best-fit 3-planet Keplerian orbital model plus one additional long-period Keplerian to model the stellar magnetic activity cycle. 
The model plotted is the one that produces the lowest while the orbital parameters annotated and listed in Tables and |4|are the 
medi an v alues of the posterior distributions, a) Full binned RV time series. Open black squares indicate pre-upgrade Keac/HIHES data 
(see ^3.1| >, open black circles are post-upgrade Keck/HIRES data, and filled green diamonds are APF data. The thin blue line is the best 
fit 3-planet plus stellar activity model. We add in quadrature the RV jitter term listed in Tablewith the measurement uncertainties for 
all RVs. b) Residuals to the best fit 3-planet plus stellar activity model, c) Binned RVs phase-folded to the ephemeris of planet b. The 
two other planets and the long-period stellar activity signal have been subtracted. The small point colors and symbols are the same as in 
panel a. For visual clarity, we also bin the velocities in 0.05 units of orbital phase (red circles). The phase-folded model for planet b is 
shown as the blue line, d) 2DKLS periodogram comparing a 2-planet plus activity model to the full 3-planet fit when planet b is included. 
Panels e) and f), and panels g) and h) are the same as panels c) and d) but for planets HD 7924 c and HD 7924 d respectively. 
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TABLE 4 

Derived Properties 


15 


^ 10 

u 
o 

O) 0 
> 

^ -5 

fC 

a: 

-10 

3000 3500 4000 4500 5000 5500 6000 6500 7000 

BJDtdb - 2450000 




BJDtdb - 2450000 



Fig. 2. — Velocity-activity correlation. Top: Binned RV time 
series of the post-upgrade Keck data with planets b, c, and d sub¬ 
tracted. Middle: Binned S'pjK time series of the post-upgrade Keck 
data only. Note the similarities between the variability in the top 
and middle panels. Bottom: Spearman rank correlation test of the 
velocities with values ( |Spearman|1904| >. 


MCM C, Te r Braak|2QQ6|) engine de scribed in Fulton et al 
(2013|) and |Knutson et al.l (|2014|). We treat the total 
KV model as a sum ot Keplerian orbits each parameter¬ 
ized by orbital period (Pi), time of inferior conjunction 
{Tcon},i)^ eccentricity (ci), argument of periastron of the 
star’s orbit (cji), and velocity semi-amplitude (Ki) where 
i is an index corresponding to each planet (b - d). An RV 
“jitter” term (crjitt) is added in quadrature with the mea¬ 
surement uncertainties at each step in the MCMC chains. 
We also fit for independent RV zero-points for the APF, 
pre-upgrade Keck, and post-upgrade Keck data. The 
long-period RV signal presumably caused by the stellar 
magnetic activity cycle is treated as an additional Ke¬ 
plerian orbit with the same free parameters as for each 
of the three planets. In order to speed convergence and 
avoid biasing parameters that must physically be finite 
and positive we step in the transformed and/or combi¬ 
nations of parameters listed in Table We add a 
penalty for large jitter values of the following form 


Anew 
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vel,n P ^jitt)’ 
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Parameter 

Value 

Units 
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0 0073 +0-0210 
U.UU/O _0.0076 


ei) sin (jJij 

-0 0031 +0-0063 

U.UUOi _0.0190 


CLb 

0.05664 

AU 

Mb sin ib 

8.68 

M0 

Sk* 

113.7 


** 

eg, 6 

825.9 lg;g 

K 

Cc COS (jUc 

0.017 


Cc sin UOc 

0.008 


CLc 

0.1134 

AU 

Me sin ic 

7.86 

M0 

Sc* 

28 35 +0-92 

^0.00 _o.89 

^0 

Teq,c 

583.6 ttZ 

K 

Cd COS UJd 

0.059 tToli 


Cd sin Ud 
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n 1 551 +0-0018 
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O.TT _o.78 

M0 

Sd* 

-|7 +0.49 
ID.i/' _o.48 

^0 

rp ** 

7 eq,d 

499 ±4 
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Stellar irradiance received at the planet relative 
to the Earth. 

Assuming a bond albedo of 0.32; t he mean 
total albedo of super-Earth size planets (|Demory| 

|2QT4t . - 


where a^e\,n is the velocity uncertai nty for each of the n 
measurements ([Johnson et al.||2011). 

We fit for the 23 free parameters by running 46 chains 
in parallel continuo usly checking fo r conve rgence using 
the prescription of Eastman et al. (2013). When the 


numb er of independent draws ['1 z as defined by Ford 
2QQ 6|) is greater than 1000 and the Gelman-Ru bin statis- 
tic (jCelman et al. 2003 Holman et al. 2006|) is within 
1% of unity lor all free step parameters we halt the fit¬ 
ting process and compile the results in Table We also 
list additional derived properties of the system in Table 
[^that depend on the stellar properties listed in Table 


4.3. False Alarm Assessment 

We attempted to empirically determine the probabil¬ 
ity that Gaussian random noise in the data could con¬ 
spire to produce an apparent periodic signal with similar 
significance to the periodogram peaks corresponding to 
each of the planets. We calculated 1000 2DKLS peri- 
odograms, each time scrambling the velocities in a ran¬ 
dom order drawn from a uniform distribution. We lo¬ 
cated and measured the height of the global maxima of 
each periodogram and compare the distribution of these 
maxima to the periodogram peak heights at the periods 
of the three planets in the original periodograms. The 
power for each periodogram within the set of 1000 is a 
Ay^ between a 2-planet plus activity model to a 3-planet 
model assuming the other two planets as “known”. Fig¬ 
ure shows that the distribution of maxima from the pe¬ 
riodograms of the scrambled data are clearly separated 
from the original peaks. None of the trials produce peri¬ 
odogram peaks anywhere near the heights of the original 
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orbital periods shorter than one day are unco mmon in 
the galaxy (aroun d 0.83 ± 0.18 % of K dwarfs; iSanchis- 


Ojeda et al.|2014), eleven ultra-short period planets have 


been foun(j^^|due to their high detectability. To search as 
carefully as possible for short period signals, we use the 
unbinned data sets, which consist of 797 RVs from Keck 
and 109 RVs from the APE. Our use of the unbinned data 
in this section also provides independent confirmation of 
the results obtained with the binned data above. 


Dawson & Eabrycky (2010) outline a rigorous proce- 
dure to distinguish between physical and alias periods. 
Our method for finding the orbital periods and distin¬ 
guishing aliases is as follows: 


1. Determine the window function of the data to un¬ 
derstand which aliases are likely to appear. 

2. Compute the periodogram of the data, determining 
the power and phase at each input frequency. 

3. If there is a strong peak in the periodogram, fit an 
N-planet Keplerian (starting with N=l), using the 
periodogram peak as the trial period. 

4. Subtract the N-planet Keplerian from the data. 

5. Compute the periodogram of the Nth planet in the 
model and compare it to the periodogram of the 
Nth planet in the data minus the model of the other 
planets. 

6. If a second peak in the periodogram has similar 
height to the tallest peak and is located at an alias 
period, repeat steps 3-5 using that trial period. 

7. If you explored an alias period, choose the model 

that minimizes best reproduces the ob¬ 

served periodogram. Subtract this model from the 
RVs. 


8. Treat the residuals as the new data set and go back 
to step 2. Examine the residuals from the N-planet 
fit for additional planets, and continue until there 
are no more signals in the periodogram. 


Fig. 3. — Top: Distribution of maximum periodogram peak 
heights for 1000 2DKLS periodograms of scrambled RV time se¬ 
ries. For each periodogram planets c, d, and the magnetic activity 
cycle were subtracted before scrambling the data set. The vertical 
dashed blue line marks the height of the original peak for planet 
b which is clearly separated from the distribution of peaks caused 
by random fluctuations. Middle: Same as the top panel for planet 
c. Bottom: Same as the top panel for planet d. 

peaks corresponding to the three planets. Because the 
2DKLS periodogram allows eccentric solutions, we ex¬ 
plore the scrambled RVs for non-sinusoidal solutions and 
are therefore sensitive to a wide variety of false alarm 
signals. We conclude that the false alarm probabilities 
for all three planets are < 0.001. 

4.4. Searching for Period Aliases 

The high cadence of the APE data set allows us to 
explore short-period orbital solutions. Traditional obser¬ 
vations on large telescopes such as Keck often yield only 
a few nights of data per year, making it difficult to deter¬ 
mine whether a short orbital period is an alias of a longer 
period, or a true physical signal. Although planets with 


The window functions of the individual and combined 
Keck and APE RV time series are shown in Eigure 
The window function is given by 

1 ^ 

= V ^ exp(-27rwij), (3) 

i=i 

where n is the frequency in units of days“^ and tj is 
the time of the jth observation. The data sets are com¬ 
plementary: 109 RVs from APE over the past year are 
well-distributed over the months and the year, and so 
are only susceptible to the daily aliases, whereas the 797 
RVs from Keck over the last decade are distributed in 
a way that gives some power to daily aliases as well as 
longer-period aliases. The power in the combined win¬ 
dow function illustrates that we might be susceptible to 
daily aliases, and weak signals in the periodogram might 
even be susceptible to monthly or yearly aliases. 

Based on a 2014 Nov 20 querv of exoplanets.org (I Wright et al.l 
pLlllHan ^ ^ ' 
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To take the perio dog:ram of the time seri es, we use a 


version of fasper (Press & Rybicki 

1989 

) written for 

Python. We find the same peaks in tJ 

tie periodogram of 


the data and residuals at 5.4, 15.3, and 24.5 days that we 
interpret as planets. The periodograms of the data and 
periodograms of the Keplerian models associated with 
these periods are shown in Figure Again, we recover 
a fourth peak at 2570 days that we attribute to long¬ 
term stellar activity due to the correlation between the 
RVs and the 5'hk values. The 3 planet plus stellar ac¬ 
tivity model found with the alias-search method yields 
periods, eccentricities, and planet masses within 1-cr of 
the results quoted in Table for all three planets, is 
consistent within l-a, but inconsistent for planets c and 
d. However, the arguments of periastron for all three 
planets are poorly constrained due to their nearly cir¬ 
cular orbits. We find a fifth peak at 40.8 days, which 
is also prominent in the periodogram of the stellar ac¬ 
tivity and is likely the rotation period of the star. The 
40.8-day signal has a strong alias at 17.1 days and so we 
test Keplerian models at both 40.8 days and 17.1 days 
to discriminate which is the true signal and which is the 
alias (see Figure]^. We find that the periodogram of the 
model 40.8 day signal better matches the alias structure 
in the periodogram of the data, and so we prefer 40.8 
days as a candidate stellar rotation period. 



Fig. 4.— The window function of the Keck and APF RV time 
series. While the APF window function has some power at fre¬ 
quency multiples of one day, it is flat otherwise, whereas the Keck 
window function has power at low frequencies (corresponding to 
long-period aliases) and power in broader swaths around the fre¬ 
quency multiples of one day. 


4.5. Chromospheric Activity 
Although this star is relatively ina ctive with val 


ues in the literat ure between —4.89 dlsaacso n fc Fischer 
2010|) and —4.85 (|Canto Martins et al.||2Ull|), some low- 
level chromospheric activity is detectable with our high- 
precision RVs. The most obvious feature is the long pe¬ 
riod signal with a period of ?^6.6 years and an amplitude 
of ~5 ms“^ . Although this signal looks promising as 
a long-period sub-Jupiter mass planet candidate, upon 
inspection of the 5'hk time series we notice that this ac¬ 
tivity indicator is highly correlated with this long-period 
signal in the velocities (see Figure [^. We interpret this 


signal as the signature of the stellar magnetic activity 
cycle of which we have observed nearly two full cycles. 
If we subtract this long-period signal from the 5 'hk val¬ 
ues and make a periodogram of the residuals we find a 
marginally significant peak at ~41 days. This is most 
likely caused by rotational modulation of starspots, be¬ 
cause 4f^_dnysj^_jma£_^eemected rotation period (38 
days, |Isaacson & Fischer||2010|) for a star of this spectral 
type and age. 

To gain confidence in the source of these Doppler sig¬ 
nals (activity or planets) we employed a time-dependent 
RV-5 'hk decorrelation that probes evolution of the RV 
periodogram on timescales of the stellar rotation. Sim¬ 
ple linear decorrelation is limited by the natural phase 
offset between RV and activity signals. At any given 
time, the measured RV and activity both depend on the 
flux-weighted fractional coverage of the magnetically ac¬ 
tive region on the stellar disk. The measured 5 'hk peaks 
when the magnetically active regions are closest to the 
disk center, where its projected area is largest and the 
star would otherwise be brightest. However, the RV 
approaches zero at the center of the stellar disk. The 
activity-induced RV shift is maximized when the active 
region is closer to the limb, ^1/8 phase fr om center (e.g. 


Queloz et al. 2001 Aigrain et al. 2012). In addition, 
the magnetically active regions continuously evolve on 
timescales of several stellar rotation periods weakening 
the RV-5'hk correlations over long timescales. 

Motivated by these limitations, our alternative RV- 
5'hk decorrelation method is effective on timescales of 
both the long-term magnetic activity cycle and the stel¬ 
lar rotation period. All APF data are excluded from this 
analysis because activity measurements are not presently 
available. 

Starting with the same outlier-removed time-series as 
described above, we bin together any radial velocity and 
'S'hk measurements taken on the same night. The pe¬ 
riod of the long-term stellar magnetic activity cycle is 
then identified from a Lomb-Scargle (LS) periodogram 
of the activity time-series. That period serves as the ini¬ 
tial guess in a subsequent single-Keplerian fit to the RV 
time-series using rv_mp_f it. Next, planet candidates are 
identified using an iterative approach in which the num¬ 
ber of fitted planets N increases by one starting at V = 1. 
V + 1 Keplerians are fit at each step, one accounting for 
the long-term stellar magnetic activity cycle: 


1. Determine period of Vth planet candidate: 

Identify the period, P/v, of the highest peak in the 
LS periodogram of the residuals to the V-planet 
Keplerian fit to the RV time-series. 

2. Obtain best V-planet solution: Perform an 
(V + l)-planet Keplerian fit with initial guesses be¬ 
ing the N Keplerians from the previous fit plus a 
Keplerian with initial period guess Pat- 

3. Repeat: Steps 1-2 for V = V + 1. 


We wish to subsequently remove residual RV signa¬ 
tures that can be confidently attributed to magnetically 
active regions rotating around the stellar surface. Since 
the lifespan of such regions is typically a few rotation 
periods, each measured RV was decorrelated with Phk 
activity measurements within ±100 days. We chose ±100 
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Fig. 5. — Left: LS periodograms of the data associated with each planet identified, from top to bottom: planet b, planet c, and planet 
d. In each panel, signals from the other planets and stellar activity have been subtracted. The phase of the frequency associated with the 
peak is given in radians. Right: LS periodograms of the best Keplerian model for each of the planets, from top to bottom, planets b, c, 
and d. The periodogram of each Keplerian model reproduces the peak period and alias structure of the data. 
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Fig. 6. — Top: a candidate periodic signal at 40.8 or 17.1 days emerges in the periodogram of the residuals to the 3-planet plus long-term 
stellar activity model. These two periods are related to each other by the one synodic month alias. Either 40.8 days or 17.1 days could 
correspond to the rotation period of the star, although the 40.8 day period is more prominent in the periodogram of the values. To 

test both periods, we model the best-fit Keplerian at 40.83 days (center) and 17.1 days (bottom) and show their periodograms, complete 
with phase information. 
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day windows because such windows were long enough to 
allow us to detect the 41 day rotation period in each win¬ 
dow, and short enough to capture long-term variability in 
the power spectrum on time-scales of hundreds of days. 
However, before discussing our decorrelation technique, 
we investigate the degree to which our RV and 5 'hk mea¬ 
surements are correlated and the reliability of our planet 
candidate signals. 

Figure shows a time-series of LS periodograms cor¬ 
responding to all RVs measured within ±100 days. Each 
panel contains more than 100 individual columns cor¬ 
responding to each periodogram derived from all data 
within ±100 days of the times for each column. These 
periodograms are only displayed at times when a mini¬ 
mum data cadence requirement (discussed below) is met 
and are grayed out elsewhere. Individual periodograms 
have been normalized independently to facilitate inter¬ 
epoch comparison of the relative distribution of power 
across all periods. In each panel from top to bottom, 
an additional planet signal has been removed from the 
data. The long-term magnetic activity cycle is removed 
from the data in all three panels. The fitted longterm 
magnetic activity-induced component has also been re¬ 
moved in all cases. After removing planet b (top panel), 
the majority of epochs/columns show significant power 
at 15.3 days, corresponding to planet c. The few excep¬ 
tions are epochs with the lowest cadence. The higher ca¬ 
dence epochs also have significant power near 24.5 days, 
corresponding to planet d. The 24.5-day signal is more 
apparent after planet c is also removed (middle panel). 
The presence of 15.3-day and 24.5-day peaks across all 
high-cadence epochs is consistent with the coherent RV 
signature of a planet. 

Figure also reveals significant RV power at ^30- 
50-day periods which can be attributed to magnetically 
active regions (i.e. spots and plages) rotating around 
the stellar surface at the ^40.8-day stellar rotation pe¬ 
riod. The activity-induced RV shifts dominate the peri¬ 
odograms after the planets b, c and d have been removed 
(bottom panel). A number of epochs also show substan¬ 
tial power at ^ 17.1 days. This is unlikely to be the 
signature of a fourth planet as the period corresponds 
to the monthly synodic alias of the 40.8-day rotation pe¬ 
riod. Moreover, a system of planets with 15.3-day and 
17.1-day periods is unlikely to be dynamically stable. 

The direct correlation between RV and stellar activity 
on stellar rotation timescales is further highlighted in 
Figure The top panel is the S'nK-analog of Figure in 
that each column is the running periodogram of all Shk 
measurements within ± 100 days. The suite of 30-50-day 
peaks mirror those in the RV periodograms of Figure 
indicating that magnetic activity is responsible for RV 
variation on these timescales. 

The observed ^20-day range in the period of these 
activity-induced peaks might seem surprisingly large, 
given that differential stellar rotation is likely of order 
several days (e.g. the observed difference in rotation pe¬ 
riod at the Sun’s equator and 60° latitude is ^5 days). 
However, some of this variation could also be statisti¬ 
cal noise, simply an artifact of the limited sampling of 
^40-day signals over a relatively small time baseline. To 
test this hypothesis, we generated a synthetic RV time 
series of the exact same cadence as our real Keck ob¬ 
servations and attempted to recover an injected signal 


of period 40.8 days. We first injected our three-planet 
model, adding Gaussian noise with standard deviation 
2.36 ms“^, the median velocity error in the post-upgrade 
Keck dataset (see Table|^. We then superposed the 40.8- 
day sinusoid of semi-amplitude 1.46 ms“^. The chosen 
semi-amplitude is the median of semi-amplitudes of all 
Keplerian fits to activity-induced RVs as determined by 
our decorrelation algorithm (see below). We applied the 
same planet-detection algorithm described above to these 
simulated data and subtracted three Keplerians from the 
best 4-Keplerian model, leaving only the 40.8-day signal. 

The bottom plot in Figure shows the resulting 200- 
day-running-window periodogram, akin to the bottom 
plot in Figure The ^10-15 day fluctuation in the 
recovered period^f the injected 40.8-day signal is consis¬ 
tent with activity-induced features in our real 5 'hk and 
RV data. We conclude that the varying power of the 
RV signal with periods of 30-50 days observed in the real 
data is a limitation of the observing cadence and not 
intrinsic to the star. 

The strong correlation observed between RV and stel¬ 
lar activity indicates tremendous potential for decorre¬ 
lation to validate our planet candidates. We apply a 
decorrelation algorithm that iterates through each RV 
measurement v{ti) where U is the time of the ith mea¬ 
surement: 

1. Define subsets Si and V^, consisting of all 5 'hk and 
RV measurements within 100 days of time 

2. Define ns and ny as the number of datapoints in 
Si and Vi respectively. If n^ and ny > ^min then 
proceed to next step, otherwise revert to step 1 for 
i = i ± 1. 

3. Identify the period, P^, of the highest peak in the 
LS periodogram of Si. 

4. Perform a single-Keplerian fit to V with fixed pe¬ 
riod Pi. 

5. Subtract Keplerian fit from v{ti). 

6. Repeat for i = i ± 1. 

As step [^indicates, we only decorrelated an RV mea¬ 
surement with Phk activity if both the number of Phk 
and RV measurements within ±100 days, ns and ny re¬ 
spectively, exceeded some minimum number, nmin- This 
was done to avoid removal of spurious correlations, n^in 
was chosen by examining the dependence of correlation 
significance on ns and ny. More specifically, for each 
RV measurement, we computed the Pearson product- 
moment correlation coefficient, r, between RV and Phk 
periodograms from all data within ±100 days. In cases 
of ns and ny between ^10-20, r values were widely dis¬ 
tributed between 0-1, suggesting correlations were spuri¬ 
ous. Supporting this notion, the best-fitting Keplerians 
to the activity-correlated RV noise had unusually large 
semi-amplitudes, as high as 10 ms“^, even in cases where 
r > 0.5. In contrast, when ns and ny exceeded 23, r 
ranged from 0.34-0.76, with a median of 0.51 and semi¬ 
amplitudes were a more reasonable 1-3 ms“^. We there¬ 
fore adopted n^in = 23 resulting in activity decorrelation 
of 115 of 281 RV measurements. 
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The top panel of Figure shows the evolution of the 
periodogram of the entire RV time series as the activ¬ 
ity decorrelation scheme is cumulatively applied to data 
points in chronological order. That is, the columnar pe¬ 
riodogram at each given time, t, corresponds to that of 
the entire RV time series, with activity decorrelation of 
all RV measurements before and including time t. Plan¬ 
ets b and c have been removed as well as the long term 
magnetic activity. The bottom panel dispays the RV pe- 
riodograms before and after activity decorrelation, cor¬ 
responding to the leftmost and rightmost periodograms 
in the top panel respectively. After decorrelation the pe¬ 
riodogram has significant power at 24.5, 35.1, and 40.8 
days and several other more moderate peaks. However, 
after decorrelation (rightmost column), the 24.5 signal 
dominates, in support of its planetary origin, while most 
other peaks, including those at 40.8 days and 35.1 days, 
have been drastically reduced, consistent with manifesta¬ 
tions of stellar activity. The single exception is the peak 
at 17.1 days. While it is unclear why this peak does not 
vanish, we remind the reader that is the monthly syn¬ 
odic alias of the 40.8 day stellar rotation period, and a 
planet with such an orbital period would likely be dy¬ 
namically unstable with planet c (P=15.3 days). The 3- 
planet result of this analysis is consistent with the other 
analyses previously described. This method of decorre- 
lating stellar noise from RV measurements will be useful 
for distinguishing planets from stellar activity in other 
high-cadence RV planet searches. 



5. PHOTOMETRY 

Photometric observations of planet host stars can be 
used to measure stellar rotation, detect planetary tran¬ 
sits, and sense false planet detections in RV signals. We 
collected 1855 differential photometry measurements of 
HD 7924 using the T 8 0.8 m automated photometric tele¬ 
scope (APT) at Fairborn observatory in Arizona. The 
measurements were collected between 31 December, 2006 
and 1 December, 2014. The instrument uses two photo¬ 
multiplier tubes to measure flux in Stromgren b and y fil¬ 
ters simultaneously. The b and y measurements are later 
combined into a single {b - 1 - y)/2 passband to improve 
signal to noise. The telescope nods between the target 
star and several comparison stars. A more detailed de¬ 
scription of the observin g procedure and da ta reduction 
can be found in H09 and [Henry et al.|(|2013|). We noticed 
season to season offsets in the photometry of HD 7924 
with a maximum amplitude of ^ 1.5 mmag. These off¬ 
sets are well correlated with the 5'hk values during the 
same time. To better characterize the photometric vari¬ 
ability on short timescales we remove these offsets for the 
subsequent analysis. 

The 1855 photometric observations are plotted in the 
top panel of Figure pR| and summarized in Table where 
we have removed yearly offsets by dividing the observa¬ 
tions in each season by the seasonal mean. We find a 
significant periodic signal in the full photometric dataset 
with a period of 16.9 days and peak-to-peak amplitude 
of 1.9 ppt. This is similar to the 17.1 day period found 
in the RVs and most likely a signature of the stellar ro¬ 
tation. Spots may form on opposite hemispheres and 
cause photometric fluctuations at one half of the orbital 
period. When the data are broken up into the individual 
observing seasons we find that the strongest photomet¬ 


Fig. 7.— Top: Periodogram of Keck HIRES RV measurements 
± 100-days versus time (JD-2450000) with our Keplerian model of 
planet b (P=5.4 days) removed. Each colored, single-pixel column 
is a unique periodogram of all RV measurements within a 200- 
day window centered on the time indicated on the horizontal axis. 
Dark grey boxes indicate times of observation having fewer than 
24 measurements of both RV and within ± 100 days. Each 
periodogram has been normalized independently by dividing all 
periodogram powers by the maximum periodgram value. Middle: 
Same as top, with planets b and c removed. Bottom: Same as top, 
with planets b, c, and d removed. 


ric periodicity happens during the 2012-2013 season and 
has a period of 41.5 days. This season corresponds to 
the time of maximal 5'hk activity index. Starspots dur¬ 
ing this time likely dominate the rotationally induced 
signal in the RV data. 

No significant periodic variability is found in the pho¬ 
tometry at the periods of the three planets. A least- 
squares sine wave fit to the data gives semi-amplitudes 
of 0.11 ppt, 0.25 ppt, and 0.17 ppt at the orbital periods 
of planets b, c, and d respectively. Figure [TT] shows the 
photometric data phase-folded to the periods of each of 
the three planets. This gives further evidence that the 
5.4, 15.3, and 24.5 day signals are planetary in nature 
and not caused by stellar activity that would be visible 
in the photometry at those periods. 

With the APT data we are able to rule out transits 
deeper than «2 ppt for planet b and 5-8 ppt for plan¬ 
ets c and d. We binned the data in bins of width 0.005 
units of orbital phase around the time of center transit 
for each of the planets. For planet b, the lowest binned 
measurement within the uncertainty window of the tran¬ 
sit time is 2.3 ± 1.2 below the mean suggesting that tran¬ 
sits deeper than 3.5 ppt do not occur. With a mass of 8.7 
M 0 we can rule out non-grazing transits of a pure hydro- 


RV power within ± 100 days 
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4000 


Fig. 8.— Top: Same as Figure]^ but for (stellar activity) 

data. Bottom: Same as Figure [^bottom panel but for a synthetic 
RV time-series created by superposing a 1.46 ms“^, 40.8-day si¬ 
nusoid with the best 3-planet model to the real Keck RVs and 
2.36 ms“^ Gaussian noise, then fitting for and removing the best 
3-planet model (P=5.4, 15.3, 24.5 days). Time sampling matches 
the real-data. Dark grey boxes indicate times of observation having 
fewer than 24 measurements of both RV and within ± 100 

days. Epoch-to-epoch fluctuation in the recovered period of the 
injected 40.8-day signal is an artifact of time-sampling, Keplerian 
fitting noise, and injected Gaussian noise. 




Fig. 9. — Top: Running periodogram of the entire RV time-series, 
where only the subset of RV measurements that precede each time 
have been decorrelated with stellar activity. Keplerian models of 
planets b (P=5.4 days), c (P=15.3 days) and the long-term-stellar- 
magnetic-activity have been removed. Dark grey boxes indicate 
times of observation having fewer than 24 measurements of both 
RV and within ± 100 days. All periodogram powers have 

been normalized to the same color scale. Bottom: Periodogram 
of entire RV time-series before and after decorrelation with stellar 
activity. The periods of peaks attributed to planets and stellar 
activity are labeled accordingly by dashed vertical lines. 


gen planet, but transits of a denser planet would not be 
detectable in our data (see Figure 11). Our constraints 
on possible transits of planets c anod are not as strong 
because their ephemerides have not yet been subject to 
intense campaigns to search for transits. However, the 
lowest binned measurements within the transit windows 
for planets c and d are 6.23 ±1.5 ppt and 3.4 ±1.5 ppt, 
respectively, and suggest that transits deeper than 7.7 
ppt for planet c and 4.9 ppt for planet d are not present 
in the data. This rules out transits of a pure hydrogen 
planet d, but the lowest binned measurement for planet c 
is consistent with the depth of a transit caused by a pure 
hydrogen planet. However, if this were a real event we 
would expect the neighboring measurements to also be 
slightly lower than the mean since the transit duration is 
longer than the width of the bins but this is not the case. 
Both of the neighboring measurements are higher than 
the mean flux level. We conclude that transits of planets 
with compositions dominated by hydrogen for any of the 
three planets are unlikely, but more and/or higher preci¬ 
sion observations are needed to exclude transits of rocky 
planets. 

6. SPITZER TRANSIT SEARGH 

A photometric campaign to look for transits of HD 
7924b using the Spitzer Space Telescope showed no evi¬ 
dence o f transiting planets lar ger than 1.16 R^{2 — a con¬ 
fidence, Kammer et al.|[2M4 ). However, they assumed a 
circular orbit which caused them to underestimate the 
uncertainty on the time of transit . We find an ephemeris 
that is inconsistent with that of|Kammer et al. (2014) 
by ^ 3 — cr, but nearly identical to that of H09. With 
the new ephemeris listed in Table we calculate a time 
of transit for the epoch during which the Spitz er ob- 
servations were collected of 2455867.07lo‘ii- The 


Kam¬ 


mer et al. (2014) observations would have covered 70% 


of the predicted ingress or egress times assuming a per¬ 
fectly edge-on viewing angle for the orbit. The a priori 
transit probabi l ity fo r planet b is 6.4%. Assuming the 
Kammer et al. (2014) observations are of sufficient pre- 
cision to detect any ingress or egress that would have 
occurred within their observational window, the a poste¬ 
riori transit probability is 2.0%. Figure 12 compares the 
observational window of the Spitzer observations to the 
predicted time of transit based on the ephemeris of this 
work. 


7. ARGHITEGTURE AND STABILITY 
7.1. Compaet Multi-planet Systems 


In order to compare the architecture of the HD 7924 
system with other compact multi-planet systems we com¬ 
piled a catalogue of similar multi-planet systems from 
Kepler. We restricted the Kepler multi-planet systems 
to those with exactly three currently known transiting 
planets in the systenR% a largest orbital period less than 
or equal to 30 days, and all three planets must have radii 
smaller than 4 R^. This left a total of 31 systems from 
the sample of confirmed Kepler systems dBorucki et al. 


20111 IBatalha et al.||2013| |lVlarcy et al.||20l4t |Kowe et al. 


2014). Their architectures are presented in Figure 13 

We notice nothing unusual about the architecture of 
HD 7924 when compared with the Kepler systems. Most 


As returned by a 2014 Nov 24 query of exoplanets.org 
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TABLE 5 

SUMMARY OF APT PHOTOMETRIC OBSERVATIONS 


Observing 

Season 

i^obs 

Julian Date Range 
(HJD - 2,440,000) 

Sigma 

(mag) 

Rrot^ 

(days) 

Full Amplitude^ 
(mag) 

2006-2007 

231 

14100-14158 

0.00212 

17.1 ±0.2 

0.0021 ± 0.0004 

2007-2008 

524 

14370-14523 

0.00215 

17.1 ±0.1 

0.0011 ± 0.0003 

2008-2009 

464 

14728-14867 

0.00208 

45.7 ±0.5 

0.0013 ± 0.0003 

2009-2010 

123 

15092-15222 

0.00212 



2010-2011 

140 

15459-15598 

0.00239 

18.1 ±0.1 

0.0024 ± 0.0005 

2011-2012 

125 

15823-15963 

0.00222 

16.8 ±0.1 

0.0021 ± 0.0005 

2012-2013 

109 

16185-16330 

0.00235 

41.5 ±0.8 

0.0027 ± 0.0006 

2013-2014 

100 

16555-16674 

0.00168 

25.1 ±0.3 

0.0022 ± 0.0004 

2006-2014 

1816 

14100-16674 

0.00197 

16.8922 ±0.0014 

0.00098 ± 0.00013 


^ Period of most significant peak in a periodogram analysis. 

2 Amplitude of best-fit sine function with the period fixed to Prot • 


of these multi-planet systems contain three planets with 
masses between 5-10 M 0 and semi-major axis between 
0.05 and 0.3 AU. The systems of Kepler-IM, Kepler- 
124:, Kepler-219, Kepler-372, Kepler-310, Kepler-127, 
and Kepler-342 all host one inner planet and two outer 
planets that are closer to each other than to the inner¬ 
most planet as in the HD 7924 system. The Kepler-372 
system is particularly similar with one planet orbiting at 
0.07 AU and a pair of outer planets orbiting at 0.14 and 
0.19 AU. This demonstrates that the HD 7924 system 
is not abnormal in the context of other known compact 
multi-planet systems. While this is not a direct demon¬ 
stration of dynamical stability, the fact that many other 
systems exist with very similar architectures gives strong 
empirical evidence that the HD 7924 system is a stable 
planetary configuration. 

The HD 7924 planetary system is one of very few RV- 
detecte d systems hosting t hree super Earths . Only HD 
40307 (|Ma!yor et al.||2009a|), and HD 20794 (|Pepe et al. 
' 2011 ) are known to host three planets with masses all 
below 10 M 0 . Both of these stars are unobservable from 
most northern hemisphere observatories. If we expand 
the mass limit to include systems with three or more 
planets with masses below 25 we find fou r additional 
syste ms that m atch these criteria, HD 69830 (Lovis et al. 
' 2006|), G1 5 81 (Mayor et al. 2009b|), HD lOTS O (iLovis 
et al.| 201 1|), and 61 Virginis (|Vogt et al.||2010|). I'hese 
systems are difficult to detect given their small and com¬ 
plex RV signals and HD 7924 is the first discovery of such 
a system since 2011 . 

7.2. Dynamieal Stability 

In order to check that the HD 7924 system is dynami¬ 
cally stable for many orbital periods we ran a numerical 
integration of the three plane t system using the MER¬ 
CURY code (|Chambers|[l999|). We started the simula¬ 
tion using the median orbital elements presented in Ta¬ 
ble and let it proceed 10^ years into the future. We 
assume that the system is perfectly coplanar with zero 
mutual inclination and we assume that the masses of the 
planets are equal to their minimum masses (sini = 0 ). 
No close passages between any of the three planets were 
found to occur (<1 Hill radius) during the entire simula¬ 
tion, suggesting that the system is stable in the current 
configuration. 

Eollowing the arguments of Eabrycky et al. (2014|) we 
calculated the sum of separations between the planets (b 


to c, plus c to d) in units of their Hill radii. Eabrycky 


et al. (2014) found that the sum of the separations tor 
the vast majority of Kepler multi-planet systems is larger 
than 18 Hill radii. The sum of the separations between 
the planets of the HD 7924 system is 36 Hill radii, pro¬ 
viding further evidence that this system is stable and not 
abnormally compact relative to the many Kepler multi¬ 
planet systems. 

8. APE VS. KECK 

Since this is the first publication from our group utiliz¬ 
ing APE data we perform a comparison between the two 
data sets to assess their relative performance. Eirst, we 
use the automated planet detection algorithm described 
in §4.1| on each data set independently. When the APE 
and Keck data are analyzed independently the 5.4, 15.3, 
and 24.5 day signals are the first to be identified in both 
cases (albeit with much lower significance in the APE 
data alone) making it extremely unlikely that these sig¬ 
nals could be caused by instrumental systematic noise. 

The binned measurement uncertainties for the APE 
data are generally higher than those of Keck (before stel¬ 
lar jitter is accounted for). This is not surprising since 
the APE data is usually collected at lower signal to noise 
than the Keck data. However, since we are allocated 
many more nights on APE we can collect more measure¬ 
ments over a shorter amount of time. It took ~13 years to 
collect 281 binned measurements with Keck, but only 1.5 
years to collect 80 binned measurements with the APE. 
In addition, HD 7924’s circumpolar track at Lick Ob¬ 
servatory and the long nights during the winter months 
in Northern California both increase the observability of 
this target. The standard deviation of the residuals to 
the best fit model are similar in the two cases at 2.8 
m s“^ for APE and 2.5 m s“^ for the post-upgrade Keck 
data. In order to compare the tradeoffs between higher 
cadence and higher precision we define a metric similar 
to the “photometric noise rate” as used by |Eulton et al.| 
( 2011 ) to compare the relative performance between pho- 
tometric data sets. We define a “velocity noise rate” as 


VNR= 

Vf’ 


(4) 


where a is the RMS of the velocities and T is the mean 
number of observations per year. We assume that no 
further signals are present in the data and/or they con¬ 
tribute negligibly to the RMS. The VNR gives the ap¬ 
proximate K in ms“^ that would be detectable with 
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Fig. 10. — Differential photometry of HD 7924 from APT. Top: 
Relative flux time-series of HD 7924 in the combined Stromgren 
(6 + y)/2 passband. Seasonal offsets are removed by dividing by 
the mean within each season. The standard deviation of the pho¬ 
tometric time series is 2.3 parts per thousand (ppt). Middle: LS 
periodogram of the photometric time series. Bottom: Differential 
photometry from APT phase-folded to the 16.9 day peak found in 
the LS periodogram and binned to widths of 0.1 phase units (red 
circles). The light grey points and open red circles show the same 
data wrapped by one period. The peak-to-peak amplitude of the 
16.9 day photometric signal is 1.1 ppt. 


S/N=l after one year of observing at typical cadence. 
Of course, we generally require S/N ^ 1 in order to con¬ 
sider a signal a viable planet candidate, but the VNR 
still gives a good reference point for comparison. We 
list the VNR and other performance characteristics for 
Keck and APF data for HD 7924 and two well-known 
RV standard stars (HD 10700, and HD 9407) in Table 
The velocities for the three stars are compared side- 
by-side in Figure [T3 The VNR for APF data is 25-50% 
smaller than that of Keck indicating that we will be sen¬ 
sitive to smaller planets once we have observed the stars 
for a comparable amount of time or, in other words, our 
APF data will be of comparable sensitivity to the exist¬ 
ing Keck RVs in 1/2 to 3/4 of the length of time it took 
to reach that sensitivity at Keck. 


We are actively tuning and perfecting our Doppler 
analysis pipeline for APF. Recently, we implemented an 
experimental technique to correct the velocities for any 
correlations of the RVs with environmental parameters 
(e.g. atmospheric pressure, CCD temperature). This 
technique reduces the RV RMS for Keck data slightly 
(5-10%), but significantly reduces the RMS for APF ve¬ 
locities in most cases and up to a factor of two in some 
cases. The RV RMS for HD 9407 and HD 10700 reduce 
to 2.24 ms“^ and 1.67 ms“^ respectively when correla¬ 
tions with non-astrophysical variables are removed. 


9. DISCUSSION & SUMMARY 

We present the discovery of two additional super-Earth 
mass planets orbiting the bright K0.5 dwarf HD 7924. 
These planets join a previously-known planet in a sys¬ 
tem with at least three super-Earth mass planets. The 
two new planets have minimum masses of Me sin 4 = 7.8 
M 0 and M^sini^ = 6.3 M 0 , and orbit HD 7924 with 
semi-major axis of Uc = 0.113 AU and = 0.155 AU. 
Both planets receive far too much radiation fro m HD 
7924 to be within th e habitable zone as defined by |Kop-| 
parapu et al. (2013) with incident stellar irradiation val¬ 
ues 114, 28, and 15 times that received from the Sun by 
Earth for planets b, c, and d respectively. Assuming that 
these planets have bond albebos similar t o the mean to- 
tal albedos of super-Earths {At = 0.32, Demory 2014) 
their equilibrium temperatures are 826, 584, and 499 K. 

The stellar magnetic activity cycle is clearly visible in 
our long-baseline RV time series and we observe nearly 
two complete cycles. We simultaneously model the RV 
shift due to the magnetic cycle and the three planets in 
order to extract accurate planetary parameters. A ten¬ 
tative RV signal from rotationally-modulated starspots 
is also found and we perform a rigorous analysis to de¬ 
termine that the planetary signals are distinct from the 
stellar activity signals. 

With the largely expanded data set we are able to re¬ 
fine the ephemeris of planet b and show that additional 
transit monitoring is needed. Since HD 7924 is near 
the ecliptic pole, it will be near the continuo us view¬ 
ing zone fo r the James Webb Space Telescope (Gardner 
et al.|2QQ6 ) and a transiting planet would make an excel¬ 
lent candidate for space-based transmission s pectroscopy. 
The T ransiting Exoplanet Survey Telescope (jRicker et al. 
2014) will observe HD 7924 once launched and will be 
able to conclusively determine if any of the three planets 
are transiting. 

This system is a good example of a compact system 
of short-period planets for which high-cadence observa¬ 
tions are incredibly valuable to determine the coherence 
of signals and detect the true physical periods as opposed 
to aliases. Since we kno w that short-period s uper-Earth 
size planets are common ( Howard et al.|2QlQ ), we expect 
our continued nearly nightly observations with the APE 
to uncover many more systems like HD 7924 from which 
we will build a comprehensive census of the small planets 
in our local solar neighborhood. 


We are very grateful for the donations of the Levy fam¬ 
ily that helped facilitate the construction of the Levy 
spectrograph on APE. Without their support the APE 
would not be contributing to the discovery of planets like 
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Fig. 11. — Differential photometry from APT of HD 7924 phase-folded to the ephemerides of the three planets. Left: Differential 
photometry phase-folded to the orbital period of planet b (top), c (middle), and d (bottom). Red circles are the photometry measurements 
grouped in bins of width 0.04 units of orbital phase. The grey shaded region spans the uncertainty in mid-transit time. The semi-amplitude 
of a sine wave least squares fit to the data is annotated at the bottom of each panel. Right: Same as the left panels zoomed in around 
the phase of center transit. The individual measurement s are omitted and the r ed circles correspond to binned photometry data in bins of 
_width 0.005 units of orbital phase. The three curves are jMandel & Agol| (|2002|) transit models for planets with the masses listed in Table 
l^and densities of iron (black), water (blue), and hydrogen fgreenj. 


TABLE 6 

APE vs. Keck Radial Velocity Precision 


Instrument 

RMS 

Median Uncertainty^ 

^obs^ Mean Cadence 

VNR3 


(ms“^ ) 

(ms“^ ) 

(days) 



HD 7924 


Pre-upgrade Keck/HIRES 
Post-upgrade Keck/HIRES 
APE/Levy 

1.20 

2.50 

2.80 

2.64 

2.36 

2.86 

7 

281 

80 

110.4 

13.5 

5.9 

0.66 

0.48 

0.36 

HD 10700 (t Ceti) 

Pre-upgrade Keck/HIRES 

2.87 

2.84 

84 

17.5 

0.63 

Post-upgrade Keck/HIRES 

2.32 

2.22 

190 

19.4 

0.53 

APE/Levy 

2.16 

2.12 

66 

7.0 

0.30 

HD 9407 

Pre-upgrade Keck/HIRES 

1.89 

1.94 

12 

162.8 

1.26 

Post-upgrade Keck/HIRES 

1.89 

1.86 

202 

18.3 

0.52 

APE/Levy 

2.34 

2.30 

104 

4.5 

0.26 


^ Stellar jitter has been added in quadrature with the binned velocities such that the of the velocities 
with respect to their median value is 1.0. 

^ Binned in units of 0.5 days. 

^ See Equation!^ 
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BJDtdb-2455000 


Fig. 12. — Comparison of th e transit search time w indow covered 
by the Spitzer observations of|Kammer et al.| (|2014j) to the transit 
times as calculated from the updated epJiemeris m this work. The 
thick black line shows the posterior probability density function 
fP DF) of the mid - transi t time for the transit that was targeted 
by|Kammer et al.| (|2014|). The blue shaded region shows the time 
window covered by their observations, and the red dashed lines 
show the most likely ingress and egress times. 

these. In their honor we informally name the HD 7924 
system “The Levy Planetary System”. We thank the 
many observers who contributed to the measurements 
reported here. We thank Kyle Lanclos, Matt Radovan, 
Will Deich and the rest of the UCO Lick staff for their 
invaluable help shepherding, planning, and executing ob¬ 
servations, in addition to writing the low-level software 
that made the automated APF observations possible. 
We thank Debra Fischer, Jason Wright, and John John¬ 
son for their many nights of observing that contributed 
to the Keck data presented in this work. We grate¬ 
fully acknowledge the efforts and dedication of the Keck 
Observatory staff, especially Scott Dahm, Greg Dopp- 
man, Hien Tran, and Grant Hill for support of HIRES 
and Greg Wirth for support of remote observing. We 
are grateful to the time assignment committees of the 
University of HawaiJ, the University of Galifornia, and 
NASA for their generous allocations of observing time. 
Without their long-term commitment to RV monitoring, 
these planets would likely remain unknown. We acknowl¬ 
edge R. Paul Butler and S. S. Vogt for many years of 
contributing to the data presented here. A. W. H. ac¬ 
knowledges NASA grant NNX12AJ23G. L.M.W. grate¬ 
fully acknowledges support from Ken and Gloria Levy. 
This material is based upon work supported by the Na¬ 
tional Science Foundation Graduate Research Fellowship 
under Grant No. 2014184874. Any opinion, findings, 
and conclusions or recommendations expressed in this 
material are those of the authors and do not necessar¬ 
ily reflect the views of the National Science Founda¬ 
tion. G.W. H. acknowledges support from NASA, NSF, 
Tennessee State University, and the State of Tennessee 
through its Genters of Excellence program. This research 
made use of the Exoplanet Orbit Database and the Exo¬ 
planet Data Explorer at exoplanets.org. This work made 
use of the SIMBAD database (operated at GDS, Stras¬ 
bourg, Erance), NASA’s Astrophysics Data System Bib¬ 
liographic Services, and the NASA Star and Exoplanet 
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Fig. 13.— The architecture of Kepler multi-planet systems com¬ 
pared with the architecture of the HD 7924 system. We plot 
the Kepler systems with exactly three currently known planets, 
a largest orbital period of 30 days or less, and a maximum plan¬ 
etary radius of 4 R^. The HD 7924 planets are plotted as filled 
black circles. The systems are sorted by the semi-major axis of 
the shortest period planet in the system. The size of the circles is 
proportional to planet radius. The radii of the HD 7924 planets 

have been calculated as = 0.371 

I-1 ® 
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Database (NStED). Einally, the authors wish to extend 
special thanks to those of Hawaidan ancestry on whose 
sacred mountain of Maunakea we are privileged to be 
guests. Without their generous hospitality, the Keck ob¬ 
servations presented herein would not have been possi¬ 
ble. 
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Fig. 14.— A comparison between RV performance at Keck and APF for two well known RV standard stars. These stars are some of the 
most well-observed stars at both Keck and APF. In the first 1.5 years of APF operations we have already collected nearly half of the data 
that has been collected at Keck over the last ~15 years. Stellar jitter has been added in quadrature with the binned velocities such that 
the of fho velocities with respect to their median value is 1.0. 
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